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We study interaction effects on the topological crystalline insulators protected by time-reversal 
(T) and reflection symmetry ( R ) in two and three spatial dimensions. From the stability analysis 
of the edge states with bosonization, we find that the classification of the two-dimensional SPT 
phases protected by Z 2 x [U(1) x T] symmetry is reduced from Z to Z4 by interactions, where 
the Z 2 symmetry denotes the reflection whose mirror plane is the two-dimensional plane itself. By 
extending the approach recently proposed by Isobe and Fu, we show that the classification of the 
three-dimensional SPT phases (i.e., topological crystalline insulators) protected by R x [U(1) x T} 
symmetry is reduced from Z to Zg by interactions. 


I. INTRODUCTION 

Recently, topological structures of gapped quantum 
states have attracted much attention. The topologically 
nontrivial phases are many-body states with gapped ex¬ 
citation spectra in the bulk and characterized by non¬ 
trivial topological structure of wave functions. A re¬ 
markable property of topological phases is the existence 
of stable gapless boundary modes, which is a source of 
various exotic properties; the boundary modes are the 
origin of quantization of the Hall conductivity in inte¬ 
ger quantum Hall systems^ and topological magnetoelec¬ 
tric effects in three-dimensional topological insulators— 
One of the important questions in this field is classi¬ 
fication of topological phases, i.e., to count how many 
topologically distinct phases exist under given symme¬ 
try. The first answer to this question is obtained for 
free-fermion systems, for which classification is summa¬ 
rized in the so-called periodic table of topological in¬ 
sulators and superconductors— Free-fermion systems 
are categorized into the ten Altland-Zirnbauer symme¬ 
try classes^ in terms of time-reversal, particle-hole and 
sublattice symmetry. The periodic table tells us that, 
in every spatial dimension, three and two out of the 
ten Altland-Zirnbauer symmetry classes have topological 
insulators/superconductors characterized by topological 
indices from Z and Z 2 , respectively. The periodic table 
of topological insulators provides helpful information for 
searching for topological materials. 

The notion of topological insulators and supercon¬ 
ductors is extended to systems of interacting fermions 
or bosons. The symmetry-protected topological (SPT) 
phases^ are many-body ground states with short-range 
entanglement which have gapped excitations in the bulk 
and gapless excitations on the boundary which are stable 
against any symmetry-preserving perturbation. The Hal¬ 
dane gap phase of integer-spin chains^ and bosonic inte¬ 
ger quantum Hall state s 10 ' 11 are examples of bosonic SPT 
phases in one and two dimensions. As for fermionic SPT 
phases, topologically nontrivial phases are also expected 
in several 4/- or 5 d- electron systems, and the Kondo 
insulator SmBg is a candidate of a topological insulator 


of correlated electrons—However, it was shown by 
Fidkowski and Kitaev that interactions can change the 
topological classification (periodic table) of free-fermion 
systems.— They demonstrated that inter-chain interac¬ 
tions in eight Kitaev chains gap out all Majorana zero- 
modes without symmetry breaking. This implies that 
topological classification of one-dimensional topological 
insulators/superconductors in class BDI is changed from 
Z to Zg by interactions. Such collapse of topological 
classification is also observed in two dimensions^— and 
three dimensions— ~ 2 - Several theoretical frameworks for 
classification of SPT phases in interacting systems have 
been developed. For example, theory based on group 
cohomolog y 7 ' 8 ' 25 or cobordisms^r— classifies possible 
topological actions for the bulk states, while theoretical 
approach using the Chern-Simons theor y 29 ' 30 or nonlin¬ 
ear sigma models^ -24 ' 31 examine the stability of gapless 
boundary modes. 

In many cases SPT phases are protected by local sym¬ 
metry such as time reversal and other internal symme¬ 
tries. However, spatial symmetry (e.g., reflection, rota¬ 
tion etc.) can also protect topological phases, as pointed 
out by Fu22 and realized in the topological crystalline 
insulator SnTe . 33 ' 34 This compound respects the time- 
reversal and reflection symmetry and has four Dirac 
cones on the (0, 0,1) surface— Having a trivial strong 
Z 2 index, SnTe is a trivial band insulator as the four 
gapless Dirac cones are not protected from opening of 
a gap by time-reversal symmetry. However, the reflec¬ 
tion symmetry about the (1,1,0) plane prohibits any 
Dirac mass term that could gap out Dirac cones, and 
the stability of surface Dirac cones is guaranteed by mir¬ 
ror Chern numbers defined on mirror planes in the Bril- 
louin zone; the topological classification is thus Z—— 
Besides SnTe, there are other candidate materials the¬ 
oretically proposed as topological crystalline insulators 
in strongly correlated electron systems: a heavy-fermion 
compound^ YbBi 2 and a d-electron system^ SrgPbO. 

Recently, effects of interactions on topological crys¬ 
talline insulators have been addressed by Isobe and Fu— 
Motivated by first-principles calculation s 41 ' 42 predicting 
that thin films of SnTe become two-dimensional topolog- 
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ical insulators (i.e., quantum spin Hall insulators), they 
studied stability of gapless edge states with internal Z 2 
symmetry that comes from reflection symmetry about 
the two-dimensional plane; they obtained Z 4 classifica¬ 
tion for thin films. Furthermore, they have developed a 
theoretical approach to classify three-dimensional topo¬ 
logical crystalline insulators which utilizes the stability 
analysis of the gapless edge states. This approach led to 
Zs classification for three-dimensional topological crys¬ 
talline insulators. Strictly speaking, however, the sym¬ 
metry of the models studied in Ref.|40is U(l)xZ 2 , where 
the U(l) symmetry denotes charge conservation^. In 
Ref. 1401 the time-reversal symmetry of topological crys¬ 
talline insulators is taken into account simply by assum¬ 
ing nonchiral gapless edge (domain-wall) modes. While 
this procedure gives correct classification, in our opinion, 
their treatment of time-reversal symmetry is not fully 
satisfying. 

In this paper we study interaction effects on topolog¬ 
ical crystalline insulators protected by time-reversal and 
reflection symmetries. That is, we classify fermionic SPT 
phases under Z 2 x [U(l) xi T ] symmetry in two dimen¬ 
sions and R x [U(l) xi T] symmetry in three dimensions, 
where T and R denote time-reversal and reflection. Using 
the models studied by Isobe and Fu, we carefully derive 
transformation laws of fields under symmetry transfor¬ 
mations and examine the stability of gapless edge modes 
under symmetry-preserving perturbations and interac¬ 
tions. We obtain Z 4 classification for two dimensions 
and Zs classification for three dimensions, in agreement 
with Isobe and Fu. 

The rest of this paper is organized as follows. In Sec. HU 
we consider a two-dimensional model respecting the re¬ 
flection symmetry whose reflection plane is parallel to the 
two-dimensional plane. We elucidate collapse of topolog¬ 
ical classification from Z to Z 4 due to interactions. In 
Sec. IIII1 by extending the argument of Ref. |40| to time- 
reversal invariant systems, we show collapse of topologi¬ 
cal classification from Z to Zs in three-dimensional topo¬ 
logical crystalline insulators. This means that eight Dirac 
cones on the surface of a topological crystalline insulator 
can be gapped out. In addition, in Sec. IIVI we point out 
that two Dirac cones can be gapped out without sym¬ 
metry breaking by attaching a fractionalized quantum 
spin Hall insulator with topological order. In Sec. [V] we 
summarize our results. 


II. TWO-DIMENSIONAL TOPOLOGICAL 
CRYSTALLINE INSULATOR 


Recent first-principles calculations predict that quan¬ 
tum spin Hall states can be realized in the (111) thin 
films of the SnTe class of three-dimensional topological 
crystalline insulators i 41 ' 42 The surface Dirac fermions on 
the top and bottom surfaces of the ( 111 ) thin films are 
gapped by intersurface coupling, and turn into a topo¬ 
logical state. Following Ref. |4j], we take the effective 


Hamiltonian for the Dirac fermions on the top and bot¬ 
tom surfaces of a thin film with an odd number of layers, 


H 2 d = J d 2 xijj\x) [( id x cr y — id y a x ) ® t z + mT x \ tp{x), 

( 1 ) 

where x — (x,y), ■j/b = (^^, V4t> ^ 24 ,)> an< 4 t ^ ie ve “ 
locity is set equal to unity. Here ip^ is the creation 
operator of Dirac fermions with spin <7 4 . on the top 

(a = 1) or the bottom (a = 2) surface. The Pauli ma¬ 
trices a 1 and r* (* = x, y, z) act on the spin and surface 
indices, respectively. The coupling of the top and bottom 
surface Dirac fermions gives the mass term mr x . 

Notice that the Hamiltonian (ID respects the time- 
reversal (T) symmetry and the local Z 2 symmetry ( g ) 
arising from the reflection symmetry with respect to the 
two-dimensional plane. Transformations of ip under the 
symmetry operations are described by 

T = - ia y K , g = ia z T x , (2) 

where K is the operator for complex conjugation. 

The presence of the Z 2 symmetry promotes the topo¬ 
logical index characterizing the two-dimensional Hamil¬ 
tonian from Z 2 (Z 2 index of time-reversal invariant in¬ 
sulators) to Z (a mirror Chern number). To under¬ 
stand this, we first note that the Chern number is de¬ 
fined for each eigenspace of g. The operators T and g 
commute with each other, [— ia v K, ia z r x ] = 0, and the 
time-reversal symmetry interchanges the eigenspaces of 
g. For example, if |+) and |—) are eigenstates of g with 
the eigenvalues +* and —i, respectively, then T|±) are 
eigenstates with the eigenvalues =F i, respectively. Thus, 
the Hamiltonian is block diagonalized into eigenspaces 
of g which are related by the time-reversal symmetry 
and characterized by the Chern numbers with oppo¬ 
site signs^. Having the integer topological number, N 0 
copies of H 2 d can have N 0 pairs of helical edge modes 
protected by time-reversal and reflection symmetries. 

To describe a Kramers pair of helical edge modes of 
the two-dimensional topological insulator, we introduce 
smooth spatial modulation in the mass term in Eq. ©■ 
This generates gapless states moving along a line where 
the sign of the mass changes. Let us replace mT x by 
m{x)T x , where m[x) takes a positive (negative) value for 
x > 0 (x < 0). This yields helical edge states localized 
at the kink (x = 0); solving the Dirac equation for the 
zeromode, 


[id x a v <g) t z + m(x)T x ] |±y) = 0, (3) 


we obtain a Kramers pair of gapless modes \+y) and | —y) 
propagating to the +y and — y directions, respectively: 


(x\±y) = exp 


r 

±ik y y — / dx m{x ) 
Jo 


l±Z/)o, 


(4a) 


with 


l±2/)o = 



1 

—i 


. (4b) 
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where |®) is the eigenstate of a;, and k y denotes the mo¬ 
mentum along the y direction. These states are trans¬ 
formed as 



In order to examine the stability of the gapless he¬ 
lical edge modes in the presence of interactions, we 
bosonize the fermionic gapless edge modes j 29 i 30 i 45 i 46 The 
Lagrandian for the bosonic fields representing the helical 
edge modes is given by 

r cIt 

£ = J — [. K^jd t (l)i{x)d x ^j(x) - d x cj) I (x)d x (j) I (x )\, 

( 6 a) 

with K = p z , where p 1 (i = x , y , z) are the Pauli matrices, 
and summation over the repeated indices I, J is assumed 
(this is assumed throughout this paper). The bosonic 
fields, 0 = ( 0 i, 0 2 ) = ( 0 + , 0 _ ), describe the edge modes 
propagating to the +y and — y directions, respectively. 
More specifically, the vertex operators : e 1 ^ : and : e 1 ^ : 
create fermions in the \+y) and | —y) states, respectively. 
The fields 0/ are defined modulo 2i r. The commutation 
relations of these fields are given by 

[0/(x),0j(y)] = inp} jSgn(x - y) + insgn{I - J). ( 6 b) 

Here, sgn(:r) equals +1, 0, and —1 for x > 0, x = 0, 
and x < 0 respectively. The second term on the right 
hand side of Eq. (l 6 bl) accounts for the anticommutation 
relation of fermions from different edge modes. Sym¬ 
metry transformations of the bosonic fields 0 / can be 
deduced from Eq. ©, which should be compared with 
Q (Q = T,g,ug). We thus obtain the follow¬ 

ing transformation rules for the bosonic fields: 

T(j){x)T~ 1 = —a x 4>(x) + 7 re 2 , (7a) 

g0(a;)g _1 = 0(a;) + |(ei - e 2 ), (7b) 

u e 4>(x)ug 1 = 4>(x) + 0(e! + e 2 ), (7c) 

where e/ is the unit vector whose Jth entry is one and 
the other entries are zero. The operators T, g , and ug de¬ 
note the time-reversal, the local reflection, and the charge 
U(l) rotation, respectively. The minus sign in Eq. (JT^,) 
is due to the antiunitary nature of T. The symmetry 
group of the system is denoted as Z 2 x [U(l) x T] in the 
notation of Ref. [29l 

We will examine whether interactions can gap out 
N 0 copies of the gapless helical modes defined above 
without symmetry breaking for N 0 = 1, • - - ,4. The 
Lagrangian for the 27Vo-component bosonic field 0 = 
(0i,..., 4>2N 0 ) T , which describes the N 0 copies of heli¬ 
cal edge modes, is given by Eq. (l 6 al) with the K matrix 


K = p z <g> 1 at 0 , 


( 8 a) 


propagating to the +y (—y) direction for odd I (even /), 
respectively. The commutation relations of the bosonic 
fields 0/ (J = 1,..., 2iVo) are given by 

[0/(x), </>j(y)\ = in(K~ 1 ) I j sgn(x - y) + wrsgn(I - J). 

( 8 b) 

The transformation rules in Eq. (ED hold if we replace 
ei and e 2 with the following 2iVo-dimensional vectors: 
ei = (1, 0,1, 0, ■ • • , 1,0) and e 2 = (0,1,0,1, • • • , 0, 1 ). 

The N 0 pairs of helical edge modes are gapped out if 
there exist potential 


/lint ^ [ Ca I dx . COS {let ' 0 T U Q ) . (9) 

Oi ^ 

with Nq linearly independent vectors l a satisfying the 
Haldane’s null vector condition^ 


l T a K~Hp = 0, (a,p = !,•■• ,2Vo) (10) 


so that the fields satisfy [l a - 0 (a:), lp-<fi{y)\ — 0 up to 27rm 
(n S Z) for a, /3 = 1, • • ■ , N 0 . The coupling constants C a 
and the phases a a are real numbers. As indicated by 
the colons in Eq. ©, the vertex operators are normal- 
ordered, 

• J2l < j hh ■ g®h01 . . gi^202 ..... g^2Ar 0 02AT O . 

( 1 !) 

where the phase factor exp(^ Yli<j is in accor¬ 
dance with the second term in the commutator in Eq. 
(l 8 bP . It makes : e~ ll cl> : the Hermitian conjugate of: e a 

There are cases when the symmetry is broken spon¬ 
taneously, even when C ln t in Eq. © is invariant under 
the symmetry transformations. The occurrence of spon¬ 
taneous symmetry breaking can be judged by finding el¬ 
ementary bosonic variables as follows. The Nq linearly 
independent vectors {li, • • ■ , In 0 } form a A^-dimensional 
lattice as R = j\l-\ + j 2 i 2 + ■ ■ ■ + jN 0 lN 0 , where j a are 
arbitrary integers (a = 1,..., No). From the integer vec¬ 
tors in the lattice, we find a set of linearly independent 
vectors {ii, I 2 , • • • ,In 0 }> from which we define a primitive 
lattice vectors {«i, • • • , uat 0 }, 


v a 


1 

gcd(L,i,--- 



( 12 ) 


where ged denotes the greatest common divisor of the 
integers in the parentheses. The vectors {«!,••• ,v No } 
form a primitive cell of the smallest volume (see examples 
below). Finally, elementary bosonic variables^ are given 
by v a ■ 0 (a = 1, • • • , No). The set of elementary bosonic 
variables {ui • 0 , • • • , vn 0 ■ 0 } are pinned to constant val¬ 
ues, when the fields {Zi • 0 , • ■ ■ , In 0 • 0 } are pinned by the 
potentail in £ int . We see that the gapless edge modes can 
be gapped out without symmetry breaking if and only if 
the set {«! • 0 , • • • , V]y 0 ■ 0 } is invariant, 

G{v i • 0, • • • , v No ■ 0 }£ -1 = {iq • 0 , • • • , v No ■ 0 }, (13) 


where ljv 0 is the Nq x Nq unit matrix. The /th compo- modulo 2 - 7 t under the symmetry transformations in Eq. 
nent of the bosonic field 07 corresponds to the edge mode ED, where Q = T,g,ug. 
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Here, two comments on the symmetry transforma¬ 
tions of the cosine terms are in order. First, the vector 
Iq = AT _1 ( 1,1, • ■ • , 1) T = (1, —1, • • • ,1, —1) is in the lat¬ 
tice generated by the elementary vectors {iq, • • • , tqv 0 }, 
when the cosine terms with {Zi, • • • , In 0 } in Ant respect 
the symmetry and gap out all the edge modes. In other 
words, the field Iq ■ 4> is pinned when all the edge modes 
are gapped out. This is understood by noting that 
:cos(2Zo •</>): is invariant under the symmetry transfor¬ 
mations. Second, it is important to take into account ad¬ 
ditional phase shifts coming from the Klein factors [the 
second term in Eq. ([30] when we examine the invariance 
of the cosine terms in A n t, in particular, under the time- 
reversal transformation, which changes the direction of 
propagation of edge modes. 

We show below that four copies of the helical edge 
modes (16al) can be gapped out without symmetry break¬ 
ing. We begin with the case of N 0 = 1. A gapping 
potential respecting the symmetry is given by 


Ant = C J dx :cos \2v ■ cf>(x)\ : 

(14a) 

with 


^ = (1,-1) T - 

(14b) 

This potential is invariant under the operations of T, g , 
and ug. Indeed, the invariance under g and ug can be 
seen as 

g (2v ■ cf>) g 1 = 2v ■ cf> + 2n, 
ug ( 2v ■ (f>) Ug 1 =2v ■ 4>. 

(15a) 

(15b) 

To verify the invariance under T, we note from Eqs. ([TJi) 
and (flTT) that 

yi . g2 iv-4> . Ji—l _ rp . g2i0i . . g — 2 i 02 . ji—1 

— ■ g2*02 . . g-2i0i . 

_ . g-2iu-<A . 

(16a) 

and vice versa, and hence 


T :cos(2u • 4>) : T _1 =:cos(2v •</>):. 

(16b) 


It turns out, however, that the ground state breaks the 
symmetry as the elementary bosonic variable v ■ <fi is 
transformed as 

g v ■ (fig -1 = v ■ <fr + 7t. (17) 

The vector v in Eq. (114b ) is the elementary bosonic vari¬ 
able respecting the U(l) symmetry. Hence, it is impossi¬ 
ble to gap out edge modes by any symmetry-preserving 
potential without spontaneous symmetry breaking when 
N 0 = 1. 

For Nq = 2, an example of pinning potentials allowed 
by the symmetry is given by 

Ant = C ^ / dx :cos[2u a • 4>(x)\: (18a) 

a—1,2 


with 

v i = (1, 0 | 0, —1), u 2 = (0,11 -1,0), (18b) 

where C is an arbitrary real number, and the vertical 
lines in tq 2 are inserted between the copies of helical edge 
modes. The invariance under g and ug can be checked 
directly as in the No = 1 case. The time-reversal trans¬ 
formation of <fi = (4> i, fo, (f> 3 , </>4) t , 

T </>T -1 = -a x © a x 4> + ir(e 2 + e 4 ), (19) 

yields 

T :cos(2tq ■ </>): T -1 = :cos(2u 2 •</>):, (20a) 
T : cos(2u 2 • cf >): T _1 = : cos(2iq ■</>):, (20b) 

where we have used Eq. CD- However, the spontaneous 
symmetry breaking occurs once the elementary bosonic 
variables {tq • <fi, v 2 ■ 4>} are pinned by A n t, as these 
variables are not invariant under g [in a similar way to 
Eq. (fl7l) ]. We note that this observation holds for ar¬ 
bitrary Z’s of the pinning potentials respecting the U(l) 
symmetry. For example, if the cosine terms with 

h = (1, —111, —1) T j (21a) 

Z 2 = (1,1|-1,-1) t , (21b) 

are used to gap out the edge modes, then the elementary 
bosonic fields {iq • 4>,v 2 ■ cf)} are also pinned (because 
li +1 2 = 2rq and l\ — l 2 = — 2v 2 ), which gives rise to 
spontaneous symmetry breaking. Thus, we cannot gap 
out edge modes without symmetry breaking. 

Let us move on to the case of Nq = 3. When all the 
helical edge modes are gapped out, the field 

Z 0 • </> = </q - </>2 + </>3 - </>4 + </>5 - <t>6, (22) 

necessarily takes a classical value (i.e., (Zo • ) = const.). 

This implies that the ground state spontaneously breaks 
the symmetry, as Iq • 4> is n °t invariant under g , 

gl 0 ■ cj)^ 1 = l Q ■ (f> + 3tt. (23) 

For No = 4, we can gap out all the edge modes while 
respecting the symmetry by taking the following cosine 
terms: 

4 f 

£| 4 t ' = C ^ / dx :cos [ v a ■ cf>(x)\: (24a) 

a=l ^ 

with 


v i = (1,0|1,0|0,—10, —1) T , 

(24b) 

^2 = ( 0 , 1 1 0 , 1 1 — 1, 0 | — 1,0) T , 

(24c) 

^3 = (1, -1 1 - 1, 1 10,0 1 0,0) r , 

(24d) 

v 4 = (0,0 I 0,0 11,-11-1, if. 

(24e) 
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The time-reversal invariance can be seen as follows: 

T :cos(t>i ■ 4 >): T -1 = :cos(i >2 ■</>):, (25a) 

T :cos(t )2 • <p) : T _1 = :cos(t»i •</>):, (25b) 

T : cos(t >3 • 4 >) : = : cos(i >3 ■</>):, (25c) 

T :cos(v 4 • <p): T _1 = :cos(t >4 •</>):. (25d) 

Invariance under the transformations of g and ug can be 
seen by straightforward calculations; 

g v a ■ (jig- 1 = v a ■ 4> (mod 27r), (26a) 

ugv a ■ (pUg 1 = v a ■ 4> : (26b) 

for a = 1, • ■ ■ ,4. When elementary bosonic variables 
v a • (jt [a = 1, ■ • ■ ,4) are pinned at some constant values 
such that (v\ ■ (p) = (v 2 ■ <p), the edge modes are com¬ 
pletely gapped out without symmetry breaking. Hence, 
we conclude that two-dimensional SPT phases form a 
Z 4 group, which is reduced from Z of the noninteracting 
fermions. 


III. THREE-DIMENSIONAL TOPOLOGICAL 
INSULATOR 

In this section, we discuss three-dimensional topologi¬ 
cal insulators with the time-reversal and reflection sym¬ 
metry. In Ref. HO, Isobe and Fu studied interaction ef¬ 
fects on topological crystalline insulators characterized 
by mirror Chern numbers. They introduced a spatially 
varying mass term to the surface Dirac Hamiltonian and 
examined the stability of one-dimensional gapless modes 
propagating along domain walls where the sign of the 
Dirac mass changes. We note that the topological crys¬ 
talline insulating phase in SnTe is an SPT phase pro¬ 
tected by both reflection and time-reversal symmetries. 
In Ref. |4(J, however, the time-reversal symmetry is not 
explicitly considered and is used only in the assumption 
of the nonchiral structure of boundary (domain wall) 
states4£ Here, we improve on their approach by keep¬ 
ing time-reversal symmetry intact at every step of calcu¬ 
lations and show that topological classification collapse 
from Z to Z§. 

Let us start with discussion on the topological invari¬ 
ant in the noninteracting case. As is well known, topo¬ 
logical phases of class All in three dimensions form a Z 2 
group. In the presence of additional reflection symmetry, 
the group structure of the topological phases is promoted 
to Zi2i When the reflection plane is the y^-plane, the re¬ 
flection operator is given by 

R = ia x P , (27) 

where o x acts on the electron spin, and P changes the 
sign of x coordinate, P : (x,y,z) —> (—x,y,z). The rel¬ 
evant symmetry group is denoted as R x [U(l) x T] in 
the notation of Ref. |29|. The topological invariant for 
insulating phases under this symmetry group is the mir¬ 
ror Chern number um = (n + — n~)/2 € Z, where nr 1 


are the Chern numbers defined in the eigenspaces of R 
on the two-dimensional fixed plane under the reflection 
in the momentum space. The time-reversal symmetry 
T = —ia y K is not closed in each eigenspace of R on the 
fixed plane. 

Following Ref.|40|, we classify the correlated topological 
crystalline insulators by introducing a spatially varying 
Dirac mass term to the surface Hamiltonian. We will 
show below that the gapless modes from eight copies of 
Dirac cones can be gapped out without symmetry break¬ 
ing. The classification of three-dimensional topological 
crystalline insulators is reduced from Z to Zs by interac¬ 
tions. 

The Hamiltonian of a single Dirac cone on the surface 
of a topological crystalline insulator is given by 

H s urf = J d 2 xxp\x){id x a y - id y a x )xp(x), (28) 


where i/d = (xpt,xp'^). Without reflection symmetry, two 
copies of the surface Dirac cones are gapped out without 
breaking the time-reversal symmetry by introducing a 
mass term in the Hamiltonian 

i?surf = (c^xxj)^ {x)[{id x a v -id y a x )®T Q +moa z ®T v ]xj){x). 

(29) 

Here r acts on the pseudospin space in which the eigen¬ 
values ±1 of t z distingsuish the two Dirac cones, and 
x /d is a four component vector, -i/d = (xjj[^, xj)\xp^, V' 24 ,); 
where xp^ a is the creation operator of a Dirac fermion 
with spin a =f, 4 . and pesudospin a = 1,2. The instabil¬ 
ity of two Dirac cones is consistent with the Z 2 classifica¬ 
tion of symmetry class AIL However, when the reflection 
symmetry is imposed, the mass term is not allowed by 
the symmetry ( Rmocr z ® t v R~ x = —mo<j z <g> t v ), and 
the two Dirac cones remain gapless; the Z classification 
follows. 

One possible way to respect both time-reversal and 
reflection symmetries while introducing a mass to Dirac 
fermions is to modulate the sign of the mass in space. 
This gives rise to one-dimensional gapless modes at the 
domain wall where the mass changes its sign. To see this, 
let us consider the following simple model: 

Hgurf = J d 2 x xp\x) [(id x a y - id y a x ) <g> r 0 

+m(x)a z ® t v }xp(x), (30) 

where m(x) takes a positive (negative) value for x > 0 
(x < 0), respectively. A sketch of the model is shown 
in Fig. [T] When the mass changes its sign as drawn in 
Fig. [jib), a Kramers pair of gapless states propagating 
along the y direction are present at the domain wall at 
x = 0 [see Fig-UJc)]. Their wave functions are written as 


(x\ ±y) = exp 


±ik y y 



dx'm{x') 


±y)o, 


(31a) 
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FIG. 1. (a) Sketch of a spatially varying Dirac mass intro¬ 

duced to surface Dirac fermions of a three-dimensional topo¬ 
logical crystalline insulator. The green dashed square denotes 
the reflection plane. We introduce a positive (negative) mass 
in the region where the red (violet) thin rectangle is attached, 
(b) Spatial dependence of the mass in Eq. (1301) . (c) Top view 
of the panel (a). Helical edge modes are localized at the kink 
x = 0. 

with 

l+y)o = *^ ® (J) , (31b) 

I -y) o = i ® , (31c) 

where | +y) and | —y) propagate to the +y and —y direc¬ 
tions, respectively. Again, \x) denotes the eigenstate of 

x. 

Under the operations of the time-reversal and reflec¬ 
tion symmetry, these states are transformed as follows: 


and 


T\+y) 


R\+y) 


= 1 - y )> 

T\-y) = -\+y), 

(32a) 

= *l +y), 

R\~y) = -i\ -y)- 

(32b) 


IV. INSTABILITY FROM COUPLING TO 
TOPOLOGICAL ORDERED PHASES 

So far, we have discussed the possibility of gapping out 
boundary modes of two- and three-dimensonal topologi¬ 
cal crystalline insulators without symmetry breaking by 
interactions among themselves. In this section, we con¬ 
sider a different situation in which gapless edge modes 
of a two-dimensional topological crystalline insulator are 
interacting with edge modes of a fractionalized quantum 
spin Hall insulator (QSHI*) In this case we find that 
even a single Kramers pair of edge modes in Eq. 0 can 
be gapped out, without spontaneous breaking of time- 
reversal and reflection symmetry. This result implies 
that two copies of Dirac cones on the surface of a three- 
dimensional topological crystalline insulator can also be 
gapped out without symmetry breaking by the coupling 
to edge modes of a QSHI*. In fact, gapped edge states 
between a QSHI and a QSHI* are discussed by Lu and 
Lee in Ref. |49|. Here we show that the same symmetric 
gapped states are obtained in the presence of additional 
reflection symmetry. 

The QSHI* is obtained by gauging the symmetry of 
the fermion number parity from a QSHI; see Appendix „ 
The Lagrangian of edge modes of a QSHI* is given by 

r /7 t * 

C = J — (-2p x I}J d t (l)id x (j)j - Vijd^rd^j) , (33a) 

with 


0 = (<^,«f) T , (33b) 


where V is a symmetric and positive definite matrix. The 
bosonic fields <p c and <j> s describe chargeon and spinon, 
respectively. These bosonic fields are transformed as fol¬ 
lows: 



Hence, introducing the bosonic fields ^ for the gapless 
modes propagating along the ±y directions, we obtain 
the transformation laws of (jr^ which are identical to Eqs. 
0 with <f> = {(j> + ,(j> ) T . Thus, proceeding in the same 
way as in the previous section, we see that four copies 
of gapless modes in Eq. (1311) can be gapped out without 
symmetry breaking by introducing cosine potentials sim¬ 
ilar to Eq. (1M1) that pin the gapless bosonic fields. Since 
the helical edge modes in Eq. (1311) are obtained by intro¬ 
ducing a domain wall in the Dirac mass that couples two 
Dirac cones, we conclude that the SPT phases of three- 
dimensional topological crystalline insulators protected 
by time-reversal and reflection form a Zs group, which is 
reduced from the Z group of the noninteracting fermions. 


Thus, the total Lagrangian of the helical edge modes at 
the interface of a two-dimensional topological insulator 
(with a mirror Chern number) and a QSHI* is found 
from Eqs. (liral) and (1351) to be given by 

[‘ (Igr 

C= — [K I}J dt(j)id x (j)j - Vj jd x (t)id x (j)j ] , (35a) 

with 

K= (~ 2 p x )ep z , V' = V@(vl 2 ), (35b) 

and 

0 = (^ S A C A + ,<I ) ~) T , (35c) 
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where I2 is the two-dimensional identity matrix, and v 
is the velocity. We seek cosine potentials of the bosonic 
fields :cos(Z • 4> + cq ): that can gap out the edge modes. 
The cosine potentials should satisfy the Haldane’s null 
vector condition in Eq. (fTOl) . I11 addition, their Z-vectors 
should be of the form Kl' with V G Z 4 , since the cosine 
potentials are composed of bosonic or fermionic vertex 
operators. The following potential terms satisfy these 
conditions: 

Qnt = C jdx [:cos(20 c - cj) + - 4>~): 

- :cos(20 s - c/) + + cf >~):] , (36) 

where C is a negative constant. These terms respect all 
the symmetry we impose: time-reversal, reflection, and 
U(l) charge. Invariance under g or ug operation can 
be seen by straightforward calculations, and the time- 
reversal invariance can be seen as follows: 

T -W) : y-i _ . e d?<t> c -<f> + -<t >-). (37a) 

j 1 - e i( 2 0 s -0 + +>Z ,_ ) : j 1 - 1 = :e ~i{' 2 4> s -4> + +<t>~) : (37b) 

In contrast to the Nq = 2 case discussed in Sec. El spon¬ 
taneous symmetry breaking does not occur in the present 
case. The elementary vectors associated with the Z vec¬ 
tors of the cosine potentials in Eq. (13(11) are given by 


= (—1, —1,1,0) T , (38a) 

V 2 = (1, —1,0,1) T . (38b) 

Applying the operator T transforms the vertex operators 
as 

T :e ivi < l , :f- 1 = - -.e™ 2 '*:, (39a) 

f-.e™* *: f- 1 = - :e ivi (39b) 

which implies that time-reversal symmetry is preserved 
as long as the pinned fields satisfy the relation 

(«1 • <t>) ~ (V2 ■ 0) = 7T. (40) 


This is indeed the case with the cosine potentials with 
C < 0 in Eq. (1361) . The invariance of the elementary 
bosonic variables Vi ■ 4> under g or ug can also be verified 
easily. Thus, we conclude that the helical edge modes be¬ 
tween a QSHI and a QSHI* are gapped without symme¬ 
try breaking even in the presence of reflection symmetry. 

In Sec. IIIII we have considered a Kramers pair of gapless 
mode induced at a domain wall of a spatially varying 
Dirac mass that couples two Dirac cones on the surface 
of a three-dimensional topological crystalline insulator. 
As a corollary of the above result, we conclude that a 
gapless helical domain-wall mode is gapped out without 
symmetry breaking when coupled to a helical edge mode 
of a QSHI*. 

V. CONCLUSION 

Motivated by theoretical proposals of topological crys¬ 
talline insulators in correlated electron systems, we 


have studied fermionc SPT phases which respect the 
time-reversal and reflection symmetry by employing the 
Chern-Simons approach. For two-dimensional systems 
respecting the reflection symmetry whose reflection plane 
is parallel to the two-dimensional plane, our analysis 
elucidates that the SPT phases form Z4 group, while 
topological classification in the noninteracting case is Z. 
Furthermore, we have addressed classification of three- 
dimensional topological crystalline insulators by extend¬ 
ing the argument proposed in Ref. Ho| with keeping time- 
reversal symmetry intact. Our analysis revealed that 
eight Dirac cones on the surface are completely gapped 
out without symmetry breaking by two-particle backscat- 
tering of gapless modes on the domain wall of varying 
Dirac mass. This leads to collapse of topological classi¬ 
fication from Z to hg. Finally, we have pointed out the 
instability of gapless boundary modes through coupling 
to topological ordered phases. 
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Appendix: GAUGING FERMION PARITY 
SYMMETRY TO FERMIONIC SPT PHASES 

In this appendix we review the fractionalized quan¬ 
tum spin Hall insulator, QSHI*, which is obtained by 
gauging symmetry of the fermion number parity Pf of a 
QSH L 49 ! 50 The effective action of the QSHI is written as 

/ / fftvp \ 

d2x y-^- K ip a i d ^ a i - - (A- 1 ) 

with K = —p z . Here, is the anti-symmetric tensor 
(e 012 = 1 , n,v,p = 0,1,2), d := {d tl d x ,d y ) 1 x := (x,y). 
The internal Chern-Simons gauge fields a/ p (/ = 1,2) 
describe the bulk low-energy excitations in the QSHI. 
Summation over repeating indices is assumed. The in¬ 
teger vector I £ Z 2 characterizes quasiparticles, whose 
equation of motion is given by 

hr = —K LJ d v a J p . (A.2) 

Equation (IA.1I) predicts gapless modes at the boundary, 
which is set to be at y = 0. The Lagrangian of the 
boundary modes is given by 

f dx 

^edge = / ^jdt^ldx^J (A.3) 

with <fi := (4^, 4>^) T [I, J = 1, 2 and (t,4-) = (1] 2)]. Here 
(j^ and r are the bosonic fields describing gapless modes 
of up-spin and down-spin quasiparticles, vf denotes the 






Fermi velocity. The commutation relations of these fields 
are given by 

[<f>i(x),^j(y)\ = -inp z i jSgn(x-y)+insgn(I-J). (A.4) 


where fields a p = Mjja p are the internal Chern-Simons 
gauge fields of this gauged system. The I\ matrix for the 
gauged system gives the mutual statistics of quasiparti¬ 
cles with integer vectors (n, n' gZ 2 ) as 


The transformation laws of these fields under T, g, ug are 
given by Eq. CD with (</>+,<(> ) replaced by Let 

us gauge the symmetry of fermion number parity. The 
corresponding operator Pf transforms the bosonic fields 
as 

PrfPj 1 = <t> + 84>p f , (A.5a) 

where 

SfpPf = -7T . (A.5b) 

Gauging this symmetry is carried out by extending the 
action as follows: 

(I) A vortex (symmetry flux) is attached to quasipar¬ 
ticles such that a quasiparticle with ( £ Z 2 going 
around a vortex acquires a phase shift l-S<j>p f which 
is equal to the phase shift from the symmetry trans¬ 
formation generated by Pf. This makes the l vec¬ 
tors of quasiparticles non-integer-valued. 


2nn T K g 1 rt! = 2irn T M T K 1 Mn'. (A.9) 

Without loss of generality, we can redefine the fields as 
a p = Xpja'J with X € GL{ 2,Z). Choosing 



we obtain the following action of the QSHI* phase: 
r / f P v P \ 

C g = J d 2 x 1^—2ffa'fatf - J . (A.ll) 

This action predicts gapless edge modes at the boundary 
y = 0, which are described by the Lagrangian 

£g,e dge = j ^ ( 2 P/, jdtfad X <!>J ~ Vgj J d x (j) I d x (j)J^ , 

(A.12a) 

with 


(II) The action (IA.1I) is extended so that the quasipar¬ 
ticles with non-integer l vectors become elementary 
excitations. 

The symmetry flux introduced in the step (I) is described 
by the vector 

lr f = (ji) ■ (A.6) 

We can see this by calculating from Eq. (IA.1I) the phase 
acquired by a quasiparticle with an l vector going around 
the symmetry flux lp f . The phase is 2nl T K~ 1 lp f , which 
should be compared with l ■ 5<j>p f . With the flux contri¬ 
bution lp f , quasiparticles are effectively described by the 
following vectors: 

l' := Mn, (A.7a) 

with 

M := ( ji n&I 2 . (A.7b) 

In the step (II), we modify the action (1A.1I) so that all the 
quasiparticles are described by integer vectors, which we 
identify with n in Eq. (IA.7I) . The modified action then 
reads 

r e pvp 

£g = J <Xx—K gi J a I ll d v d J p - (A.8a) 

with 

^ = /=i ’ 2 ’ ( A - sb ) 


0 T = («/> s ,-r). (A.12b) 

Here, V g is a positive-definite symmetric matrix. This 
edge Lagrangian is equivalent to the Lagrangian in 
Eq. (l33l) . The bosonic fields (4> s ,(j) c ) are related to 
(0 f , <P l ) as follows: 

( r, -r) ++ ( <t>\ &) m{x t )-\ (a. 13 ) 

or equivalently, 

£)«($;£)• ,A ' 14) 

From this correspondence, we obtain the transformation 
law of fields <p s and <j> c in Eq. (1511) . It follows from the 
correspondence that 

T :e 2i ^“: T~ l (A.15a) 

f : e 2irf,c : T _1 -O- — : e i ^ + ^ ) : . (A. 15b) 

One can see how the time-reversal operator transforms 
the vertex operators in the pinning potentials in Eq. (1361) 
in the following way. With the correspondence, 

:eiW +_ 0 -_ 2 ^). ^ ;e i( 0 +-<r^W ):i (A.16a) 
;e i(^++ 0 --^ . e i(b + +^^ t -0 i ) ;) (A.16b) 

applying the time-reversal operator yields 

f ;e Tb + -<r- 2 <n. f -1 ^ :e iW + -b“-^ t +b J ') ;) (A.17a) 
T • e*(^ ++ W ~ 2( t >c ) • f- 1 ^ +(/>~-<l> f 17b) 
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from which Eq. (l36l) follows. Here the bosonic fields (f> + 
and 4>~ describe the helical modes of a QSHI. In a sim¬ 
ilar way, we obtain Eq. (13911 by noting that the vertex 


operators have following correspondence: 

(A.18a) 

■.e™ 2 *: ++ :e i (* _ ~* 4 ‘):= -i (A.18b) 

which implies the following relations: 

T :e™ 1 ^: f -1 O - ::, (A.19a) 

f : : f- 1 o - : : . (A.19b) 
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